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The variability of sea surface conditions has been observed from 
a low flying aircraft by a laser wave profiling system for fetch 
limited wind speeds of 25 M/SEC in the North Sea. Wave profiles 
obtained with the laser system have been analyzed and show that wave 
growth occurs simultaneously at all frequencies and that an equilibrium 
value for the higher frequency components is eventually reached, but 
not before substantially higher (overshoot) values are obtained. 
Simultaneous photography of the surface has been analyzed and show 
that 32% of the surface is covered with white caps, foam and streaks. 
This result is in good agreement with a semi-empirical relationship 
incorporating both the wind speed and the local wave spectrum which 
predicts 26% white water for the conditions observed. 
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INTRODUCTION 


For a number of years various scientific organizations have been 
conce r nod with developing techniques to remotely observe ocean surface 
parameters of wave height and wind speed in a quantative manner. The 
purpose of these efforts is to provide information about the "state- 
of-the-sea" in support of operational requirements and to improve 
environmental forecasts. 

During March of 1969, a study of the active and passive micro- 
wave characteristics of the ocean surface under the influence of 
various wind speeds and wave heights was staged out of Shannon, 

Ireland. This experiment was a cooperative effort of several scienti- 
fic organizations and included three aircraft and scientists from the 
Naval Oceanographic Office, NASA Goddard Space Flight Center, NASA 
Manned Spacecraft Center, NASA AMES Research Center, New York University, 
and the University of Kansas. The three aircraft were equipped with a 
variety of remote sensors capable of measuring the microwave brightness 
temperature, radar backseat ter ing cross section, surface temperature 
and wave heights. 

Between March 6 and March 13, five coordinated flights were con- 
ducted over Atlantic weather stations "I" and "J" which are located 
within about 720 KM of Shannon and are routinely occupied by British 
and French weather ships. During these flights the highest winds 
encountered averaged about 16 m/sec, an unfortunate development as 
the sites were selected due to the high probability of occurence 
of winds greater than 20 m/sec and seas greater than 8 meters. The 
weather situation was particularly interesting to the local Irish 
meteorologists who marveled at the "mildest weather in twenty years!" 

As time was running out a weak meteorological disturbance formed 
just south of Ireland and moved slowly eastward while a building high 
pressure system moved eastward from Iceland into Norway between 
March 12 and March 14. As a result of these developments, the pressure 
gradient in the North Sea gradually increased and by 2300 GMT, March 13, 
offshore winds in the North Sea reached 20 m/sec. 

By early morning, 14 March, it became apparent that these winds 
would hold, or increase, during the day and plans were made to conduct 
a flight experiment between the coast of Denmark and Scotland, the 
area of highest winds. 
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The plan of the experiment was to observe the behavior of the 
microwave signature of the surface under steady state wind conditions, 
but varying (growing) wave conditions. The NASA AMES Research Center 
Convair 990, a four engine jet aircraft, equipped with a scanning 
19.35GHz microwave radiometer, a laser wave profiler, and a medium 
resolution infrared radiometer was utilized to obtain data from a 
distance of 160 kilometers off the coast of Denmark to the vicinity 
of the Shetland Islands. The NASA MSC P3A aircraft equipped with an 
active radar scatterometer was utilized to obtain radar cross section 
data at a point midway between the coasts of Denmark and Scotland. 

This report describes a portion of the results obtained with the 
laser wave profiler aboard the Convair 990. 


INSTRUMENTATION 


The Laser Wave Profilometer 


The laser wave profilometer is a standard model Geodolite 3A 
airborne altimeter manufactured by Spectra Physics Inc., Mountain View, 
California (Figure 1). Figure (2) is a schematic diagram of the 
Instrument. 1 The ranging technique consists of amplitude modulating 
a continuous wave helium-neon laser of red light centered at 6328 A. 

The reflected light is collected by an 8 inch Schmidt Cassegrain 
telescope, detected by a photo multiplier, amplified and phase compared 
with the modulation frequency of the transmitted beam. The phase 
difference between these two signals is proportional to the transit 
time of the light, and hence the range. This range is directly 
proportional to the wave length of the modulation frequency. The 
highest modulation frequency available is 49.17MHz with a wave length 
of 3.048 meters. A total of five modulation frequencies are available 
yielding range selections of 10*-, 10^, 10-*, 10^ and 10^ feet. By 
means of a range extender circuit each of these selections may be 
extended depending on the variability of the terrain. For wave 
measurements only the full scale increments of 3.048, 6.096 and 
30.48 meters are used.* 


Static tests of the laser system have been conducted aboard an 
offshore tower. 2 Figures 3 and 4 are excerpted from reference (2) 
and show good agreement with a reference resistance wire wave staff. 

Initial flight tests of the system were conducted off Atlantic 
City, New Jersey, and in the vicinity of ARGUS Island, a U. S. Navy 
Research Tower located near Bermuda, BWI, in 60 meters of water. 


* For more details the reader is referred to reference (1). 
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Figure 5 is an example of profiles of surfaces waves as they shoal and 
break on the beach shown in conjunction with simultaneous strip 
photography. 3 Figure 6 presents a comparison of averaged wave spectra 
obtained from three two-minute tracks flown near the ARGUS Island 
Tower with that derived from a twenty-minute sample of wave measure- 
ments obtained from the tower's resistance wire wave staff. 

Since the aircraft is a moving reference, it is necessary to 
convert the observed wave spectrum to fixed coordinates. The con- 
version technique used involves accounting for the speed of the 
aircraft relative to the phase speed of each wave frequency component 
and assumes that all waves are traveling in the direction of the 
wind. ^ Since this is an approximation, the presence of swell and 
the spreading of wave energy with distance can lead to errors. In 
the case presented, little swell was present and the errors associated 
with directional spreading appear to be minimal for most purposes. 

Another source of error is that of aircraft motions of heave, 
pitctv,and roll. Figure 7 depicts the energy spectrum of heave motions, 
mapped to fixed coordinates, of the type of aircraft involved (a 
Lockheed Super Constellation) during the ARGUS Island experiment when 
the surface winds were 12 m/sec. It can be seen that essentially 
all of the aircraft heave (vertical) motions are concentrated at 
equivalent wave frequencies of less than .07HZ. The time series of 
roll angles experienced during this same three-minute track averaged 
about + 1.5°. For flight altitudes of 200 M this corresponds to 
range (wave height) errors on the order of .15 meters. From these 
data, it can be seen that the majority of range errors associated with 
aircraft motions are small within the wave frequency pass band. As 
a result, the most convenient technique for removing these errors is 
simply high pass filtering, where the low frequency cut off selected 
depends upon the turbulence and, for roost cases, is less than that 
associated with a true wave frequency of .07HZ. For extremely precise 
measurements, additional effort must be expended. For the data 
presented in this paper, the filtering technique was employed with 
apparently good results. 


METEOROLOGICAL SITUATION 


As noted in Section I, the general synoptic situation was a 
result of a weak disturbance moving eastward from Ireland simultaneous 
to an eastward displacement of a high pressure system located over 
Iceland. These two events caused a gradual tightening of the gradient 
in the North Sea and corresponding surface winds increased slowly, but 
steadily, to a maximum average of 25 m/sec by 1200 GMT, 14 March, as 


64-5 


reported by the German beacon vessel "LHHT", located at 57.7°N, 

03°E (Figure 8). 

Since only one of the reporting ships shown in Figure 8 was 
equipped with an anemometer (the rest being Beaufort estimates) it 
was felt necessary to derive the surface wind field by constructing 
isotachs of the geostrophic wind from a fine scale analysis of the 
surface pressure field. 

Surface pressure data we® obtained from a considerable number 
of British and continental shore sites and isotachs of the geostropic 
wind were drawn relying heavily on the ship reports. Near the ship 
"LHHT", a ratio of surface (20m) to geostropic winds of 70% gave 
good agreement although in nearby regions of large shear, higher geo- 
stropic wind ratios are apparent. This is reasonable considering 
the effects of lateral turbulent momentum transfer. The above analysis 
considered the effects of stability after the technique of Cardone.* 

In order to verify the validity of the inferred surface wind 
speeds, table 1 was constructed. This table shows wind speeds com- 
puted from the Litton model LTN 51 inertial system aboard the air- 
craft. The flight level winds reported are one minute averages com- 
puted approximately every five minutes during the lower altitude tracks. 
These values were averaged over twenty minutes of flight time (shown 
in brackets) and reduced to the equivalent twenty meter anemometer 
height assuming neutral stability and a logarithmic wind profile 
extending to an altitude of at least 43.3 meters. 7 The rather con- 
sistent agreement between the wind speeds derived by this technique, 
the geostropic surface winds and the ship reports, lends considerable 
confidence to the presumed wind fields. 


OBSERVATIONS OF WAVE GROWTH 


The intent of the experiment was to begin the downwind flight 
track as near to the upwind shore as possible and then to fly down- 
wind until the wave spectrum was essentially n fully developed". This 
was not entirely practical, however, as clearance procedural problems 
of the aircraft prevented approaching the shoreline within less than 
about 157 kilometers. Further, due to previous limited success be- 
cause of relatively low wind speeds, the reported winds were greeted 
with some scepticism. As a result, the flight track (dashed line of 
Figure 8) was planned to best accommodate several possible meteorological 
situations which might be present, and events seem to confirm the final 
plan selected as a good compromise. 
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Figure 9 presents plots, at their appropriate fetches, of all 
available measurements of significant wave height, Hg*. The fetch 
was measured as simply the distance from the coastline in the direction 
of the implied surface wind direction. It is seen that H<< generally 
increased with fetch to about 370KM but decreased thereafter. As a 
reference, a theoretical relation (Inoue-Cardone) is plotted for 
22 m/sec, a suitable wind speed for much of the data as determined 
from the analysis of all available wind inputs.^ The wind speed for 
the particular location of each He plotted was obtained from the 
geostropic analysis and is shown below each value. Up to 370 kilo- 
meters the variation of Hg with fetch is in good agreement with the 
theoretical relation, which applies for infinite duration, fetch 
limited seas. Beyond 370 KM, the measured wave heights are considerably 
below the theoretical relation and depart increasingly with fetch. 

This latter behavior is generally explainable by duration effects, 
as the latter measurements plotted are from a wind field that ex- 
hibited considerable variability in wind duration. In particular, 
the durations are considerably longer in the near shore (southern- 
most portion of the first downwind run) area. An inspection of the 
previous six-hourly sea level pressure analysis indicate no signifi- 
cant variation in the pressure gradient there for the 24 hours pro- 
ceeding the flight. Considering that the winds prior to 24 hours 
were only slightly lower, these durations are probably sufficient for 
the seas to be fully developed at their respective fetches. From 
the vicinity of "LHHr 1 northwestward, however, the durations are 
extremely limited. The 3 hourly wind reports at "LHHT", Figure 10, 
suggest that the duration of 20 m/sec winds was only about 12 hours 
and that duration of 23 m/sec winds (or higher) was only about 3 hours. 
Toward the end of the last downwind run, durations were probably more 
limited, considering the way in which the synoptic scale pressure 
gradient was changing. 

The behavior of the data below 370 KM suggests that these data 
might be employed to study growth, as the wind field there might be 
considered to be reasonably steady and homogeneous. 


Results 

A convenient way to study the gross behavior of a developing sea 
in this case is presented in Figure 11 which shows the plots of spectral 
density (M^-sec) as a function of frequency (true) and fetch. This 
contour analysis is based upon spectra computed for 160 KM, 176 KM, 

203 KM, 231 KM, 268 KM, and 323 KM fetches. Though this resolution 


* Hg is defined as the average value of the highest one third waves. 


is not as great, the analysis is conveniently comparable to those 
obtained in the radar altimeter fetch-limited study at lighter winds by 
Barnett and Wilkerson. The analysis in Figure 11 is considerably 
smoother than those of their study, largely due to the decreased 
resolution. The major features to be deduced from this figure are: 

1. At frequencies above about .13 Hz, the spectral density 
does not change significantly with fetch. Thus these spectral 
components are already saturated, or fully developed, and in 
equilibrium with the wind field. 

2. Below about .13 HZ, all spectral components are increasing 
in intensity with fetch simultaneously. 

3. The spectral peak generally moves toward lower frequencies 
(as indicated by dashed line) with increasing fetch, in general 
agreement with existing concepts of spectral development. A 
rather surprising feature, however, is that the magnitude of 
the spectral density at the peak does not increase markedly with 
fetch beyond about 185 KM. 

4. A given spectral component will overshoot its eventual 
equilibrium value. That is, if the development of a single 
frequency component is followed fetchwise, it is seen that the 
spectral density increases rapidly to a maximum value at a 
given fetch and then decreases gradually to a lower equilibrium 
value at greater fetch. .The effect is most noticeably between 
frequencies of .095 and .115. The overshoot effect at higher 
frequencies is not observed but presumably occurred at shorter 
fetches (160 KM) as the data suggest that the fetch at which the 
effect occurs increases with decreasing frequency. This striking 
feature of the North Sea spectra is perhaps the most remarkable, 
as it implies the operation of mechanisms in wave generation not 
yet treated satisfactorily be existing theories. The "overshoot" 
effect has been reported in wave tank studies and in the field 
experiment of Barnett and Wilkerson and its appearance here lends 
support to the fact that this is a "real" phenomena. The 
magnitude of the overshoot effect observed in the Barnett and 
Wilkerson study, that is, the ratio of spectral density at eventual 
equilibrium to that at the peak of the overshoot, was scattered 
between .35 and .75. For the North Sea mission, this ratio could 
only be computed for frequencies of .111Hz, .108Hz, and ,105Hz 
where values of .38, .50, and .41 occurred. 

Some of the features discussed above are evident if the spectral 
density is plotted versus fetch for selected frequencies. Thus 
Figure 12a shows a saturated frequency. Superimposed upon the data 
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is the spectral behavior for the Inoue-Cardone growth theory which is 
based upon a modified Miles-Phillips resonance instability growth 
mechanism and the Pierson-Moskowitz fully developed spectrum. ^ it 
is seen that generally in 12 a-f, the equilibrium value observed is 
30-40% less than that indicated by theory. This is largely explain- 
able by the different equilibrium range constant of .52X10-2 found 
for the North Sea spectra as compared to .8X10*2 used in the Pierson- 
Moskowitz spectra. The overshoot effect is clearly evident in Figure 
12 b-h. The theoretical growth, of course, does not indicate over- 
shoot, as it is based upon the relation (in the early stages of 
growth) 


d (S) - A+BS (1) 

dt 

where S is spectral density, t is time, A is the resonance mechanism 
parameterization and B is the instability mechanism parameterization. 
Dissipation effects are included by allowing the growth rate to de- 
crease as the fully developed spectral value is approached. For 
infinite duration limited fetch growth, (1) is written: 

C 0 ds(f) - A+BS (f) 
dx 

where X is the fetch, and Cg is the group velocity of each wave 
frequency, f. Disregarding the effect of the overshoot, it may be 
said that the Inoue-Cardone growth (which predicts the variation of 
Hs fairly well) is too fast for frequencies above .078Hz and too 
slow for frequencies below .069Hz. The rapid growth of these 
frequencies is surprising, especially considering that these spectral 
components should be under the influence of dissipation by bottom 
friction. 

It is possible to quantitatively compute the values of A and B 
from plots 12, a-o, provided dissipation effects can be considered 
to be small (not considering bottom friction) if the spectral density 
is less than 30% of its eventual equilibrium value*. Within this 
restriction, the computation can be performed on frequencies .087, 
.083, .078, .074, .069, and .063 Hz and will be performed on spectra 
recomputed for greater fetchwise resolution. 


* After Barnett and Wilkerson 


WHITE CAP COVERAGE 


A semi-empirical theory for the dependence of foam cover on wind 
speed and sea conditions has been preposed by Cardone.^ An important 
part of the theory is the model used to describe the growth of the 
wave spectrum due to energy transferred from a turbulent wind profile- 
this model is also used in the wave forecasting schemes developed at 
New York University. 


The growth rate of a spectral wave component S( f ) in terms of 
the wind field specified at 19.5 meter s above the sea surface, and the 
friction velocity Undefined as where t is the surface stress 


and p is the air density, may he described as 


dsiil - 

dt 


A(f,Ui9. 5 ) + B(f,U*)-S(f) [ l-f<I> 


( i-s ml 
\ S ~ 


( 2 ) 


where S(f) - spectral intensity at frequency 
A ■ Phillips Resonance growth term 
B ■ Miles-Phillips Instability Growth term 
Sqo ■ Pierson Moskowitz Fully Developed Spectrum 


The value of the term U + depends upon the local wind field and the 
air-sea temperature difference. All non-linear dissipative mechanisms 
that would act to limit growth are modeled implicitly by the use of 
the term 


1 - 


<P> 

b oo 


which causes each spectral value to approach its equilibrium or sat- 
urated value given by the Pierson-Moskowitz formula. S(f) which is a 
function of wind speed and frequency and does not take into effect the 
overshoot phenomena discussed earlier. 

The Phillips resonance term can be thought of as an initial 
excitation mechanism - necessary to start growth for parts of the 
spectrum with no energy initially - but small in comparison to the 
instability growth term which is responsible for nearly all the energy 
transfer. 

Energy transfer does not cease after a spectral component reaches 
its equilibrium value, and the energy thus transferred can be thought 
of as being dissipated in wave breaking. Whitecaps are a manifestation 
of wave breaking and thus the whitecap coverage may be closely related 
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to the energy transfer to the fully developed portion of the spectrum. 
This energy transfer is given by the expression: 

00 

E = P w • g *y*B*S*6*df 
o 

where 6 =1 when S=Soo 
6 =0 when S<Soo 


This hypothesis was tested with the whitecap data obtained 
photographically by Monahan x on the Great Lakes. The spectral growth 
model was used to hindcast the wave spectrum for each case reported by 
Monahan from corresponding estimates of stability, wind speed, and 
fetch (durations were unknown and assumed to be infinite). The energy 
transferred to the fully developed portion of the spectrum was then 
correlated with the observed whitecap percentage. The resulting linear 
regression is shown in equation 4. 

Wp ■ - .0185 + .893x10“^ .E(ergs/cm^/sec) (4) 


where Wp ■ Fresh Water foam density in percent 
E ■ Energy dissipated in breaking waves. 


Salt water whitecaps are apparently more persistent than fresh 
water whitecaps due to differences in the bubble-size spectra and hence 
for similar conditions of wave breaking, whitecap coverage should be 
expected to be greater over the oceans. An approximate salt water 
correction factor of 1.5 has been suggested by Monahan from results of 
a laboratory study of comparisons of fresh water and salt water white- 
cap data.* 


An example of the computation of percent whitecaps is shown in 
Figure 13 for a wind speed of 25 if./ sec and fetch of 180 n.mi. The 
dashed line is the fully developed spectrum for 25 m/sec, the solid 
curve is the hindcasted spectrum for a fetch of 180 n.mi. and the 
dash-dot line shows the contribution of the factor B *S to the energy 
dissipation computation. The significant wave height is calculated to 
be 26 feet. The salt water whitecap cover is computed to be 267.. 

This result compares favorably with the value of 327. calculated from 
photography obtained during the North Sea experiment. 

Figure 14 is an example of white cap conditions for the 25 m/sec 
winds observed during this flight. Note the large amount of surface 
area covered by thin streaks oriented in the direction of the wind. 
Figure 15 is an example of white-cap density for a 12 m/sec wind speed. 
The most striking difference between the two is in the amount of 
streaking evident. This characteristic of the wind driven ocean is 
important when considering the relationship of remotely obtained micro- 
wave signatures of the ocean to the surface wind and wave conditions. 
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CONCLUSIONS 


By the use of Photographic and Laser observations of sea surface 
conditions it has been found: 

1. Existing wave theory can be used with good results to predict 
the rate of growth of the significant wave height with increasing fetch. 

2. A migration of the spectral peak toward lower frequencies with 
increasing fetch is observed. 

3. All wave frequencies are receiving energy simultaneously. 

4. The spectral density of a particular frequency component will 
overshoot its eventual equilibrium value by as much as 50%. 

5. Semi-empirical calculations of white cap density based on the 
rate of energy transfer and the state of development of the wave 
spectrum agree well with observations for the nearly fully developed 
case. 


It should be noted that wave and white cap data such as discussed 
herein are very difficult to obtain. Thus extension of these results 
to other conditions of wind speed, fetch, and stability must be done 
with caution. Experimental efforts are continuing in this area and 
the emphasis will be on the short fetch and stability criterion and 
the characteristics of the surface driven by extremely high winds such 
as associated with hurricanes. These studies will be conducted both in 
a laboratory wave tank and in the field and should contribute signifi- 
cantly to future theories describing the wind driven seaway. 
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TABLE I 

A Comparison of Winds Obtained from the NASA CV990 over the North Sea and 20-Meter Level Winds 
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Fig. 4. Comparison of wave-energy spectra derived from laser and 
wave-staff recordings of surface waves. 
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WHITE WATER CREST 
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Fig. 5. Airborne laser profiles of surface waves obtained in conjunc- 
tion with simultaneous continuous strip photography. 
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. 6. Comparison of wave-power spectra obtained with the 
laser profiles aboard a U. S. Navy aircraft with that obtained 
from an in-situ wave staff. 
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Fig. 7 Aircraft heave displacement spectrum showing energy 
contribution in surface wave pass band 
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Fig. 8 Ship reports in the North Sea 1200 GMT, March 14, 1969. 

Dashed line is the flight track of the NASA CV990. The digit adjacent 
to each report is the wind speed to the nearest knot within the range 
indicated by the bar or flag. 
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(575°N, 03°E ) ON MARCH 14,1969 AT 3 HOUR INTERVALS. 
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Fig. 11. Contours of equal spectral density (in 2 -s.) 
on a frequency-fetch diagram. 
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Fig. 12. Plots of spectral density versus fetch 
for selected frequencies. 
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Fig. 13 


Pig. 14 25 m/sec vindspeed. 

Note extensive streaking present which 
accounts for 277 of the total white water 
content of 317. 



Fig. 15 12 m/sec windspeed. 

Note that thin streaks account for 
only 37. of the total of 107 white 
water. 



